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PREFACE 


Any  detailed  study  of  the  temperature  in  the  lower  atmosphere, 
as  well  as  heat  balance  studies  at  the  earth’s  surface,  must  neces¬ 
sarily  consider  conditions  at  the  surface  itself:  e.g.,  the  prop¬ 
erties  of  the  ground  regarded  as  a  medium  for  the  transfer  of  heat. 

This  report  examines  the  computation  of  ground  surface  tempera¬ 
ture  in  the  Rand  two-level  general  circulation  model  with  specific 
reference  to  the  inclusion  of  soil  heat  flux.  Comparative  results 
of  various  numerical  integrations  are  presented  to  determine  the  most 
realistic  technique  for  the  computation  of  the  ground  surface  temper¬ 
ature.  It  is  hoped  that  this  discussion  will  contribute  to  the  im¬ 
provement  of  the  accuracy  of  computed  ground  surface  temperature  in 
the  simulation  or  climatic  charges. 

Related  Rand  reports  on  the  dynamics  of  climate  are  R-877-ARPA, 

*  Documentation  of  the  Mintz-Arakaua  T  >00 -Level  Atmospheric  General 
Circulation  Model  (December  1971);  P.-1318-ARPA,  The  Radiation  and 
Heat  Budget  of  the  Mintz-Arakaua  Model:  January  (October  1973);  and 

R  - J.544-ARPA ,  The  Radiation  and  Heat  Budget  of  the  Mintz-Arakaua  Model: 
July  (forthcoming) . 
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SUMMARY 


The  computation  of  ground  surface  temperat 


- r-.-fcure  in  the  Rand  two-leve' 

atmospheric  general  circulation  model  has  been  examined  with  specific 

reference  to  the  inclusion  of  soil  heat  capacity  as  well  as  soil  heat 
flux. 

In  the  control  experiment,  a  global  ground  surface  temperature 
distribution  was  calculated  by  solving  the  heat  balance  equation  at 
the  earth-atmosphere  interface,  assuming  a  zero  soil  heat  capacity  as 
wel,  as  zero  heat  flux  into  the  soil  by  conduction.  The  results  based 
on  48-hour  integration  of  the  model  for  Janua.y  show  the  following: 

1.  The  daily  maximum  surface  temperature  is  too  high,  especially 

in  equatorial  and  tropical  regions  of  both  the  summer  and  win- 
ter  hemispheres. 

2.  The  amplitude  of  the  diurnal  temperature  oscillation  is  un¬ 
realistic  in  these  regions. 

3.  There  is  no  time  lag  between  the  maximum  solar  radiation  and 
the  maximum  ground  surface  temperature. 

The  computation  of  ground  surface  temperature  (bare  land  only)  was 
repeated  with  the  same  conditions  (e.g.,  zero  heat  capacity)  as  those 
in  the  control  experiment  eXoept  that  an  explicit  parameterized  formu¬ 
lation  for  soil  heat  flux  was  included  in  the  heat  balance  equation. 

The  results  of  these  experiments  showed  relatively  lower  amplitudes  of 
daily  maximum  surface  temperature  as  well  as  diurnal  range.  However, 
the  maxima  in  soil  heat  flux,  solar  radiation,  and  ground  temperature 
occurred  at  the  some  local  time.  This  feature  is  at  variance  with  ob¬ 
servations,  which  show  that  the  maximum  in  soil  heat  flux  precedes  the 

maximum  in  solar  radiation  by  2  hours  and  the  maximum  in  ground  tempera- 
ture  by  3  hours . 

The  model  was  then  integrated  to  compute  ground  temperature  by 
solving  a  prognostic  equation  that  included  nonzero  soil  heat  capacity 
as  well  as  soil  heat  flux.  The  results  of  this  experiment  showed  the 
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most.  realistic  distribution  of  ground  surface  temperature  with  respect 
to  amplitudes  of  ground  temperature  and  diurnal  range,  as  well  as  phase 
relationship  between  the  ground  temperature,  solar  radiation,  and  soil 
heat  flux. 
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I .  INTRODUCT ION 


SIGNIFICANCE  OF  GROUND  SURFACE  TEMPERATURE 

The  heat  flux  and  the  temperature  distribution  at  the  ground  sur¬ 
face  and  in  the  soil,  down  to  a  depth  where  temperature  waves  are 
damped  to  a  negligible  amount,  are  subjects  of  considerable  importance 
to  many  meteorological  problems.  Meteorologically,  the  ground  surface 
(as  the  lower  limit  of  the  atmosphere;  is  an  intermediary  for  the  en¬ 
ergy  exchange  between  the  atmosphere  and  the  underlying  surface.  One 
of  the  important  heat  sources  aifecting  the  atmospheric  circulation  is 
the  transfer  of  sensible  and  latent  heat  between  the  surface  and  the 
atmosphere.  The  direction  of  the  energy  flow  across  the  interface  is 
primarily  determined  by  the  difference  between  the  ground  surface  tem¬ 
perature  and  the  air  temperature  immediately  above  it.  Thus  in  studies 
of  atmospheric  circulation  it  is  necessary  to  know  the  ground  surface 
temperature  accurately. 

FACTORS  DETERMINING  THE  GROUND  TEMPERATURE 

The  observed  temperature  field  at  the  ground  (and  within  the  soil) 
is  the  result  of  a  number  of  complex  factors.  The  factors  can  be 
grouped  into  two  general  categories:  the  external  and  the  intrinsic. 
The  external  factors  are  meteorological  elements,  such  as  radiation, 
rain,  snow,  air  temperature,  humidity,  and  heat  transfer  by  the  wind. 
Although  the  relative  importance  of  each  of  these  varies,  the  radia¬ 
tion  balance  generally  is  dominant.  The  intrinsic  factors  are  the 
properties  of  the  soil  itself,  such  as  the  texture,  moisture  content, 
topography  (including  slope  and  altitude),  and  vegetation  cover.  There 
are  also  the  special  effects  of  an  urban  environment  and  of  fires.  It 
is  evident  that  to  understand  the  ground  temperatures,  both  the  exter¬ 
nal  and  intrinsic  factors  and  their  interrelationships  must  be  taken 
into  consideration. 


RESEARCH  STATUS 


Observational 

The  literature  contains  extensive  descriptions  of  the  surface 
temperature  characteristics  based  upon  observations .  These  include 
the  variations  of  the  surface  temperature,  the  amplitude  and  time  of 
maximum  ground  temperature,  and  the  effects  of  surface  cover  on  the 
surface  temperature.  Geiger  (1965)  has  discussed  these  fact-rs  in 
considerable  detail  in  his  treatise  The  Tlimate  ,7<?_  ’  the  Ground.  How 
ever,  it  is  also  seen  that  observational  records  are  quite  set rce  and 
the  distribution  of  ground  temperature  measuring  stations  over  the 
earth's  surface  xs  very  uneven.  There  is  no  universally  uniform  prac¬ 
tice  of  measurement ;  most  frequently  the  methods  of  observation  are 
determined  by  needs.  As  a  consequence,  the  representativeness  and  ac¬ 
curacy  of  available  data  are  uncertain,  and  the  results  based  on  these 
data  are  usually  approximate  and  not  necessarily  mutually  consistent 
(Sutton,  1953). 

Theoretical 

In  the  field  of  energy-transfer  climatology,  there  is  a  long  his¬ 
tory  of  modeli'g  based  upon  analytical  solutions  to  sets  of  differen¬ 
tial  equations  (see,  for  example,  Ingersoll  et  al.,  1948;  Brunt,  1934; 
Lettau,  1951;  Eckert  and  Drake,  1959;  Lonnouist,  1962). 

In  solving  the  heat  balance  equation  at  the  earth-soil  interface, 
conduction  is  the  most  important  process  of  heat  transfer  in  the  soil. 
Although  soil  is  not  a  homogeneous  medium  and  there  are  other  mecha¬ 
nisms  of  heat  transfer  in  the  soil,  the  Fourier  theory  of  heat  conduc¬ 
tion  does,  in  most  cases,  apply  to  real  conditions  in  the  atmosphere 
(Chang,  1958).  Though  the  rigorous  theoretical  treatment  of  heat 
conduction  in  soil  is  outside  the  scope  of  this  study,  we  have  given 
in  the  Appendix  a  general  outline  of  the  procedure  to  obtain  the  ana¬ 
lytic  solution  to  the  heat  conduction  equation  (Eq.  (A-3))  as  well  as 
an  expression  for  soil  heat  flux  (Eq.  (A-6)). 

The  analytic  solutions  of  the  heat-budget  equation  of  the  earth's 
surface  have  been  discussed  by  Lettau  (1951),  who  has  developed  a 
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theoretical  model  of  surface  temperature  oscillations.  This  model 
gives  the  amplitudes  and  phase  lags  of  both  the  diurnal  and  annual 
courses  in  terms  of  external  conditions  and  the  physical  properties 
of  the  soil  and  atmosphere.  The  solutions  of  this  and  various  other 
analytical  models  indicate  the  following  characteristics  of  the  sur¬ 
face  temperature. 

1.  For  zero  heat  capacity  there  is  nc  time  lag  between  the 
ground  temperature  and  the  soil  heat  flux. 

2.  For  non-zero  heat  capacity,  the  soil  heat  flux  leads  the 
surface  temperature  by  3  hours  for  diurnal  course  and  1.5 
months  for  annual  cycle. 

3.  The  amplitude  of  surface  temperature  also  depends  upon  heat 
capacity;  for  zero  heat  capacity,  the  amplitude  cannot  show 
any  appreciable  decrease. 

In  more  recent  meteorological  research,  the  ground  surface  tem¬ 
perature  is  usually  computed  by  solving  the  earth-atmosphere  interface 
heat  balance  equation  with  the  assumption  of  continuity  of  the  inter¬ 
face  temperature.  These  computations  are  generally  made  as  a  part  of 
more  complex  atmospheric  circulation  models-both  boundary  layer  models 
(Estoque,  1963;  Pandolfo  et  al. ,  1965;  Sasamori,  1970)  and  general  cir¬ 
culation  models  (Manabe  et  al. ,  1965;  Gadd  and  Keers,  1970;  Gates  et 
al. ,  1971;  Delsol  et  al.,  1971;  Kasahara  and  Washington,  1971).  The 
basic  assumptions  are  that  the  heat  capacity  of  the  earth  is  zero  and 
that  there  i'  no  net  heat  flux  through  the  air-earth  interface.  Re¬ 
cently  it  has  been  proposed  in  some  general  circulation  models  (Corby 
et  al.,  1972;  Arakawa,  1972)  to  compute  the  ground  temperature  by  solv¬ 
ing  a  prediction  equation.  Myrup  (1969)  and  Outcalt  (1972)  have  de¬ 
veloped  a  digital  surface-climate  simulator  model  to  study  diurnal 
surface  thermal  and  energy  transfer  regimes.  Their  models  are  based 
on  the  equilibrium  temperature  theory,  which  states  that  given  a  set 
of  astronomical-temporal  atmospheric  and  surface  boundary  conditions, 
there  is  one  and  only  one  surface  temperature  which  balances  the  en¬ 
ergy  conservation  equation  across  the  surface  of  the  earth.  Jacobs 
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and  Brown  (1973)  use  a  heuristically  derived  iterative  technique  to 
solve  the  quartic  form  of  the  heat  balance  equation  in  order  to  obtain 
an  equilibrium  surface  temperature. 

The  ground  temperature  values  obtained  in  the  above-mentioned 
works  have  been  discussed  at  length  in  the  literature;  the  results 
necessarily  reflect  the  assumptions  of  the  model.  For  example, 

Sasamori  (1970)  has  reported  that  the  predicted  surface  temperatures 
are  higher  than  observed;  he  attributes  this  discrepancy  to  the  as¬ 
sumption  of  no  horizontal  temperature  advection.  Estoque  (1963)  also 
reported  his  model's  tendency  to  overestimate  the  soil  heat  flux  dur¬ 
ing  most  of  the  day  but  to  underestimate  it  about  midday.  The  general 
circulation  models  (Gates  et  al. ,  1971)  indicate  highly  unrealistic 
values  of  the  ground  surface  temperature  as  well  as  the  diurnal  varia¬ 
tion,  especially  in  African  and  South  American  regions. 

Computational 

Apart  from  some  variations  in  the  formulations  of  the  components 
of  the  heat  balance  equation  among  different  models,  the  most  note¬ 
worthy  characteristics  of  the  incorporation  of  heat  flux  in  tne  soil 
are  as  follows: 

1.  While  a  number  of  local  atmospheric  and  soil  boundary  layer 
models  include  soil  heat  flux  explicitly ,  only  Delsol  et  al .  (1971) 
have  done  so  in  a  general  circulation  model. 

2.  Some  general  circulation  models  (Gates  et  al.,  1971;  Manabe 
et  al.,  1965)  have  completely  neglected  the  soil  heat  flux  as  well  as 
soil  heat  capacity. 

3.  Other  general  circulation  models  use  a  parameterized  technique 
to  compute  soil  heat  flux.  For  example,  Kasahara  and  Washington  (1971) 
specify  that  soil  heat  flux  is  three-tenths  of  the  sensible  heat  flux 
into  the  atmosphere;  Gadd  and  Keers  (1970)  consider  soil  heat  flux  to 
be  a  fraction  of  the  net  radiation  flux.  Corby  et  al.  (1972)  who  use 

a  prognostic  equation  include  soil  heat  capacity  in  the  equation  but 
do  not  compute  soil  heat  flux  explicitly. 
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PURPOSE 

purpose  of  this  research  is  to  examine  the  computation  of  ground 
surface  temperature  In  an  atmospheric  circulation  model  with  specific 
reference  to  the  inclusion  of  sol!  heat  flux.  We  propose  to  use  the 
existing  version  of  Rand's  two-level  general  circulation  model  to  compute 
the  ground  temperature.  Various  numerical  integrations  will  be  performed 
by  using  different  formulations  of  soil  heat  flux,  and  by  solving  a  time 
dependent  equation  for  the  ground  temperature  which  includes  soil  heat 
flux  as  well  as  soil  heat  capacity.  The  results  will  be  compared  to  de¬ 
termine  the  most  realistic  method  for  the  computation  of  the  ground  sur- 
face  temperatures. 
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_I_I.  THE  BASIC  CIRCULATION  MODEL 


The  basic  model  equations  used  in  this  study  are  those  described 
in  a  comprehensive  description  and  documentation  of  a  two-level  atmo¬ 
spheric  general  circulation  model  by  Gates  et  al.  (1971).  The  original 
model  of  Mintz  and  Arakawa  has  been  adapted  by  Rand  to  simulate  global 
climate,  as  well  as  to  test  atmospheric  predictability,  and  has  under¬ 
gone  a  series  of  modifications  and  improvements.  F.gure  1  shows  a 
schematic  cross  section  of  the  two-layer  model.  The  model  equations 
are  set  in  the  so-called  c  coordinate  system  in  which  the  earth's  sur¬ 
face  is  always  the  coordinate  surface  o  =  1.  Here  o  is  the  dimension¬ 
less  vertical  coordinate  given  by  *  =  (p  -  p^/^  _  p  )#  where  p  ig 

the  pressure,  pt  is  the  (constant)  tropopause  pressure,  and  p  is  the 
(variable)  surface  pressure.  The  results  of  numerical  experiments 
with  the  model  have  been  described  and  discussed  by  Gates  (1972). 

What  follows  is  based  on  these  two  Rand  reports. 
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The  model  computes  horizontal  wind  velocity,  temperature,  and 
pressure  at  each  of  the  two  levels  1  and  3,  which  approximate  the 
400-mb  and  800-mb  surfaces.  The  mixing  ratio  is  presently  computed 
.kIj  at  the  lower  level  3.  The  vertical  motion  is  given  by  5;  it  is 
determined  at  a  central  level  (  600  mb)  and  is  prescribed  to  be  zero 
at  the  top  of  the  model  atmosphere  (•  *  0,  p  200  mb)  and  at  the 
earth's  surface  (c  =  1). 

The  diabatic  heating  in  the  model  is  provided  by  the  net  radia¬ 
tion,  the  release  of  latent  heat,  and  sensible  heat  transfer  from  the 
surface.  The  radiative  heat  transfer  is  itself  affected  by  water 
vapor  and  clouds,  with  the  latter  determined  as  a  result  of  both  large- 
scale  and  convective  processes.  The  moisture  source  is  provided  by 
the  evaporation,  and  the  sink  is  the  precipitation  released  in  both 
the  large-scale  and  cumulus  scale  condensation. 

A  thin  boundary  layer  (with  subscript  4)  is  assumed  at  the  air- 
ground  interface.  This  in  effect  implies  that  the  quantities  with 
subscript  4  apply  to  the  anemometer  level  in  the  atmosphere.  The 
parameters  at  the  ground  surface  are  signified  by  the  subscript  g.  It 
may  be  remarked  here  that  the  quantities  at  level  4  (top  of  atmospheric 
boundary  layer)  are  not  predicted;  instead  they  are  obtained  from  the 
predicted  values  of  the  quantities  at  levels  1  and  3.  The  ground  sur¬ 
face  temperature  is  determined  by  solving  the  heat  balance  equation 
at  the  air-ground  interface.  Over  the  oceans,  the  sea  surface  tempera¬ 
ture  is  prescribed  and  held  fixed  throughout  the  simulation. 

The  grid  spacing  in  the  horizontal  is  4°  latitude  and  3°  longitude. 
The  finite  difference  approximations  have  been  described  in  detail  by 
Gates  et  al.  (1971).  The  time  integration  proceeds  in  steps  of  6 
minutes  (=  At);  the  forcing  terms  of  the  governing  equations  are  com¬ 
puted  every  half  hour  (i.e.,  every  fifth  time  step).  For  all  other 
relevant  details  about  the  basic  circulation  model,  the  reader  is  re¬ 
ferred  to  Gates  et  al.  (1971)  and  Gates  (1972). 


n 
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III.  SOLUTION  OF  THE  HEAT  BALANCE  EQUATION  WITHOUT 
SOIL  HEAT  FLUX;  THE  CONTROL  EXPERIMENT 


As  indicated  in  the  Introduction,  we  planned  a  series  of  experi¬ 
ments  to  determine  the  effect  of  soil  neat  flux  on  the  ground  surface 
temperature.  Accordingly,  we  first  integrated  the  basic  circulation 
model  without  incorporating  the  soil  heat  flux  in  ground  surface  tem¬ 
perature  computations.  This  is  designated  as  the  control  experiment, 
l.tus  for  the  control  experiment  we  computed  the  ground  temperature  by 
solving  the  heat  balance  equation 

R  +  H  +  LE-S  =  0  (1) 

at  the  earth-atmosphere  interface.  Here  R  is  the  net  long-wave  radia¬ 
tion  emitted  from  the  surface;  H  and  LE  are,  respectively,  the  sensi¬ 
ble  and  latent  heat  fluxes  from  the  surface;  and  S  is  the  total  solar 
radiatior  absorbed  by  the  ground.  As  described  by  Gates  et  al.  (1971), 
the  heat  flux  components  (except  for  S)  depend  upon  ground  temperature 
>,both  at  time  t  +  At  as  well  as  time  t) .  Thus,  throughout  this  study, 
we  use  backward  implicit  method  for  computing  ground  surface  tempera¬ 
ture.  In  this  study,  since  there  is  more  emphasis  on  comparison  be¬ 
tween  different  experiments  rather  than  absolute  results  of  any  one, 

.'e  have  used  exactty  the  same  formulations  for  H  and  LE  as  in  Gates 
fit  al.  (1971).  Thus  substituting  for  R,  H,  and  LE  in  Eq.  (1)  and  re¬ 
arranging,  we  obtain  the  following  expression  for  computing  the  ground 
temperature  at  time  t  +  At: 


(omitting  superscript  t  from  terms  on  the  right-hand  side).  Here 
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C 


H 


o .  C  C 
4  D  p 


V 
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where  ^  =  density  of  air  at  anemometer  level  (level  4) 

CD  =  surface  drag  coefficient 
Cp  =  specific  heat  of  air  at  constant  pressure 
V's  =  measure  of  surfate  wind  speed 

Tg  =  (aPPearin8  right-hand  side)  ground  surface  temperature 
at  time  t 

Also  =  air  temperature  at  anemometer  level 
L  ■  latent  heat  of  evaporation 

qs,q4  =  mixin8  ratios  at  the  surface  and  in  air  jusL  above  the 
ground 


GW  -  dimensionless  measure  of  ground  wetness  taking  values  from 
0  (for  dry  conditions)  to  1  (for  completely  wet  conditions) 
The  sign  convention  preceding  the  symbols  in  Eq.  (1;  is  chosen 
so  that  it  is  positive  when  the  flux  is  away  from  the  ground  and  nega¬ 
tive  when  it  is  toward  the  surface.  Equation  (1)  is  formulated  ac¬ 
cording  to  the  law  of  conservation  of  energy,  assuming  that  the  earth's 
surface  contains  no  heat  but  that  a  considerable  exchange  of  heat  oc¬ 
curs  across  it.  if  may  be  noted  that  Eq.  (1)  also  assumes  horizontal 
homogeneity  so  that  only  vertical  heat  transfer  is  considered;  hori¬ 
zontal  advection  is  of  considerable  importance  in  nature  when  two  dif¬ 
ferent  kinds  of  surfaces  exist  side  by  side  and  where  ground  is  pur.  to 
ferent  uses.  lor  the  purposes  of  our  experiment,  we  apply  Eq.  (1) 
over  bare  land,  snow-covered  land,  or  ice-covered  land.  For  ice- 
covered  ocean,  the  surface  heat  balance  equation  is  modified  by  the 
addition  of  a  term  to  include  conduction  of  heat  tnrough  the  ice;  over 
the  ocean  itself  we  prescribe  the  sea  surface  temperature. 


With  initial  and  boundary  conditions  appropriate  to  January,  the 
two-level  model  was  integrated  for  a  period  of  48  hours.  The  ground 
temperature  was  calculated  by  using  Eq .  (2).  The  results  of  the  control 
experiment  have  been  compared  with  other  experiments  by  comparing 
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daily  maximum  surface  temperatuies  and  the  diurnal  range  of  surface 
temperature  at  all  non-oceanic  points  of  the  grid.  For  more  detailed 
analysis,  we  selected  three  representative  land  points  over  the  globe, 
as  shown  in  Fig.  iy  which  also  shows  the  land  areas  of  the  globe  as 
resolved  by  the  model. 

F  .gure  3  shows  the  global  distribution  of  i daily  maximum  ground 
surface  temperature  as  simulated  by  the  model  for  the  month  of  January. 
The  realism  of  this  distribution  cannot  be  ascertained  for  want  of 
global  daily  ground  surface  temperature  maps.  However,  an  effort  will 
be  made  to  draw  inference  on  the  basis  of  available  information. 

MAGNITUDE  OF  MAXIMUM  SURFACE  TEMPERATURE 

There  are  very  few  records  of  maximum  surface  temperature  for 
natural  surfaces,  although  estimates  have  jeen  made  of  the  maxima  of 
surface  temperature  as  a  function  of  the  air  temperature  at  the  con¬ 
ventional  thermometer  shelter  level.  For  example,  Johnson  and  Davies 
(1927)  estimated  that  in  the  hottest  areas  of  the  globe  the  surface 
temperature  cannot  exceed  90°C.  This  was  based  on  an  air  temperature 
of  60°C  at  a  height  of  20  cm  above  the  ground.  However,  there  is  no 
direct  evidence  of  any  natural  surface  temperature  as  high  as  this. 
Also,  considering  that  convective  processes,  which  have  a  regulating 
influence  on  the  temperature,  become  more  efficient  for  high  surface 
temperatures,  it  is  doubtful  whether  substantially  higher  surface  tem¬ 
peratures  can  in  fact  be  reached.  These  estimates  also  do  not  take 
into  account  other  factors  that  influence  the  surface  temperature, 
such  as  the  effect  of  water  content  on  the  thermal  conductivity  of  the 
soil  and  the  movement  of  air  over  the  surface. 

In  view  of  the  above,  it  is  reasonable  to  expect  that  maximum 
"natural"  surface  temperatures  generally  do  not  exceed  60°C  in  the 
winter  hemisphere.  Consequently  we  infer  that  the  control  experiments 
results  show  an  unrealistic  distribution  of  ground  surface  temperatures 
in  the  continental  areas  south  of  30°N,  i.e.,  in  the  tropical  region 
of  winter  hemisphere  and  in  the  entire  summer  hemisphere.  Specifically, 
in  the  Sahara  region,  the  distribution  is  more  typical  of  late  spring/ 
early  summer  rather  than  January  (in  May,  the  highest  temperature  of 


tONGIIUDI 


Land  areas  as  resolved  by  the  model  grid 
locations  of  selected  points  A,  B,  and  C 


LONGITUDE 


ground  surface  temperature  (°C  )  -  Control  experiment  (January) 
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3  C  was  recorded  in  Sahara)  (Geiger,  1965).  Similarly,  temperatures 
of  90°C  in  equatorial  Africa  appear  to  be  too  high;  since  this  i  a 
very  wet  and  forested  region,  it  is  highly  unlikely  that  such  a  high 
value  for  maximum  surface  temperature  could  be  attained  in  nature. 

Physically,  the  maximum  and  minimum  surface  temperature  is  reached 
when  the  flow  of  heat  into  the  soil  exactly  balances  the  flow  outward; 
the  actual  value  attained  depends  not  only  on  the  radiation,  evapora¬ 
tion,  and  heat  transfer  xn  the  air  but  also  on  the  heat  transfer  in 
the  soil.  Since  we  have  not  incorporated  the  heat  transfer  in  the 
soil  in  Eq.  (1),  it  seems  likely  that  this  exclusior  is  largely  re¬ 
sponsible  for  the  "unrealistic"  distribution  of  ground  surface  tem¬ 
perature  found  in  the  control  experiment. 


DTURNAL  RANGE  OF  GROUND  SURFACE  TEMPERATURE 

Figure  4  shows  the  global  distribution  of  the  diurnal  range  of 
the  surface  temperature  obtained  in  the  control  experiment.  In  nature 
the  amplitude  of  the  surface  diurnal  variation  changes  considerably, 
particularly  with  the  time  of  the  year.  In  winter,  the  daily  tempera¬ 
ture  oscillation  at  the  surface  is  much  smaller  than  in  summer  (Fig. 

5).  Other  observations  (Geiger,  1965)  show  that  at  Pavlosk  (60°N)  the 
diurnal  range  in  summer  is  15°C  while  in  winter  it  is  about  2  or  3°C. 

In  Central  Europe  (Vienna)  the  amplitudes  for  winter  and  summer  are 
3  C  and  22  C,  respectively.  In  general,  outside  the  tropics,  for  clear 
days  and  bare  surface,  25°C  is  a  fairly  representative  value  of  diurnal 
oscillation  in  summer  (Sutton,  1953).  In  winter  this  value  is  about 
2  or  3°C .  It  can  be  seen  from  Fig.  4  that  the  diurnal  oscillation 
values  for  the  control  experiment  are  rather  unrealistic,  exceeding 
50  C  in  South  America.  Also,  within  summer  tropical  regions,  obser¬ 
vations  (Sinclair,  1922)  have  shown  a  diurnal  range  of  56°C  in  the 
extreme,  whereas  the  control  experiment  shows  a  diurnal  range  in  ex¬ 
cess  of  70°C  over  wide  areas  of  tropical  regions,  even  in  the  winter 
hemisphere.  Thus,  although  in  general  the  diurnal  oscillations  ob¬ 
tained  in  the  control  experiment  are  not  realistic,  the  model  at  least 
attempts  to  include  the  nature  of  the  surface  and  the  state  of  the  sky, 

both  of  which  affect  the  amplitude  of  diurnal  variation  of  the  surface 
temperature. 
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The  tine  of  maximum  surface  terrperature  is  also  of  considerable 
interest  in  analysis  of  heat  transfer  in  the  atmosphere.  Observations 
have  shown  that  on  clear  days  the  surface  attains  its  maximum  tempera¬ 
ture  about  1  hour  after  the  time  of  maximum  solar  radiation  (Lonnquist, 
1962 ;  Sutton,  1953,  p.  197).  It  may  be  mentioned  here  that  this  lag 
of  the  daily  maximum  temperature  increases  with  depth  in  the  soil. 
Figure  6  shows  a  plot  of  computed  ground  temperature  anc.  solar  radia¬ 
tion  as  a  function  of  local  time  at  two  selected  locations  (A  and  B) 
of  the  grid.  It  is  seen  that  there  is  no  time  lag  between  maximum 
solar  radiation  and  maximum  ground  temperature.  This  result  also  fol¬ 
lows  from  the  analytical  solutions  and  is  tc  be  expected  when  heat 
capacity  is  neglected. 

In  summary,  the  results  of  the  control  experiment  show  that 

1.  The  daily  maximum  surface  temperature  is  too  high,  especially 
in  equatorial  and  tropical  regions  of  both  the  summer  and 
winter  hemispheres. 

2.  The  amplitude  of  the  diurnal  temperature  oscillation  is  un¬ 
realistic  in  these  same  regions. 

3.  There  is  no  time  lag  between  the  maximum  solar  radiation  and 
the  maximum  ground  surface  temperature. 
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IV.  SOLUTION  OF  HEAT  BALANCE  EQUATION  INCLUDING 
SOIL  HEAT  FLUX:  EXPERIMENTS  1  TO  3 


In  the  previous  section  we  have  described  the  ground  surface  tem¬ 
perature  as  calculated  in  the  control  experiment  with  the  Rand  two- 
level  atmospheric  circulation  model.  In  this  section  we  describe  the 
results  of  a  series  of  numerical  experiments  which  were  performed  to 
determine  the  effect,  on  the  ground  surface  temperature,  of  incorpo¬ 
rating  the  soil  heat  flux  at  the  surface,  assuming  zero  heat  capacity 
of  the  ground  itself.  (The  situation  when  the  heat  capacity  is  not 
zero  is  considered  subsequently  In  Sec.  V.) 

In  these  experiments  the  ground  surface  temperature  was  obtained 
as  a  solution  of  the  generalized  heat  balance  equation 

R  +  H+LE+G-S  =  0  (3) 

where  R,  H,  LE,  and  S  have  the  same  meaning  and  formulations  as  in  the 
control  experiment  (Eq.  (1)),  and  G  is  the  soil  heat  flux  at  the  sur¬ 
face.  Equation  (3)  was  solved  for  bare  land  points  only.  Table  1 
shows  the  formulation  of  G" used  in  the  various  experiments. 

Experiment  1 

The  rate  at  which  heat  flows  through  soil  at  a  depth  z  is  deter¬ 
mined  by  the  thermal  conductivity  and  the  temperature  gradient  at  that 
level.  Thus  at  z  =  0  (surface),  the  soil  heat  flux  G  may  be  computed 
from  the  expression 


G  = 


where  is  the  (constant)  thermal  conductivity.  Strictly  speaking, 
this  equation  cannot  be  applied  to  real  soils  because  they  are  not 
homogeneous.  But  for  most  purposes  it  holds  when  an  appropriate  value 
of  is  chosen  (de  Vries,  1963). 
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Table  1 


FORMULATIONS  OF  SOIL  HEAT  FLUX  G  USED  FOR  COMPUTING 
'.ROUND  SURFACE  TEMPERATURE  WITH  EQ.  (3) 


Experiment 

Formula  for  Soil 
Heat  Flux  C 

Re  f  erence 

Control 

G  =  0:  no  flux 

1 

°  '  D  <TC  -  V 

Estoque  (1963),  (L)  ; 

Pandolfo  (1965),  (L)  ; 

Myrup  (1969),  (L) ; 

Delsol  et  al .  (1971)  , 

(GCM) ;  Sasamori 

0970),  (L) 

2 

| 

G  =  URN 

Gadd  and  Keers  (1970) , 

(GCM) 

3 

G  =  Q.3H 

Kasahara  and  Washington 

- _ 

LE  =  H 

0971),  (CCM) 

NOTE:  In  reference  column,  L  =  local  atmospheric 

and  soil  boundary  layer  models;  CCM  =  atmospheric  gen¬ 
eral  circulation  models. 

SYMBOLS:  ‘  =  thermal  conductivity 

=  ground  surface  temperature 

D  =  damping  depth 
l'D  =  soil  temperature  at  depth  D 

R^  =  net  radiation 

U  =  G/V  empirical  parameter 

Since  our  circulation  model  does  not  have  explicit  levels  within 
the  soil,  we  approximate  the  above  expression  for  C  as 


G  =  > 


(4) 


where  Tg  is  the  temperature  of  ground  surface,  and  TD  is  the  tempera¬ 
ture  of  soil  at  depth  D  (at  which  the  diurnal  variations  become  negli¬ 
gible!  .  The  values  of  Tg  which  can  be  obtained  by  using  Eq.  (4)  in 
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Eq.  (3)  are  dependent  upon  the  values  assigned  to  >,  ,  and  D.  Table 

2  shows  the  values  of  these  quantities  used  in  some  numerical  models. 
It  may  be  noted  that  Eq.  (4)  is  an  approximate  form  of  Eq.  (A-6)  of 
the  Appendix  if  we  neglect  the  heat  capacity. 


Table  2 

VALUES  OF  D,  TQ ,  >  USED  IN  SOME  NUMERICAL  MODELS 


Model  Type 

D  (cm) 

td  (°K) 

X  (cal  cm  1  sec  3  deg  3) 

Reference 

Local  boundary  layer 

50 

294.14 

1.44  x  io-3 

Estoque  (1963) 

50 

295 

0.8  x  io-J  to  1.9  x  io“3 

Pandolfo  (1965) 

20 

294.3 

2.6  x  io  3  (rural  areas) 

Myrup  (1969) 

_  T 

9.3  x  10  (urban  areas) 

General  circulation 

500 

280 

2.1  x  io-3 

Delsol  et  al. 

(1971) 

The  thermal  conductivity  A,  which  is  very  important  in  estimating 
the  soil  temperature,  depends  upon  the  porosity,  moisture,  and  organic 
matter  content  of  the  soil  (Chang,  1958).  De  Vries  (1963)  has  shown 
that  a  can  be  estimated  from  theoretical  considerations  as  a  function 
of  the  moisture  content  of  the  soil.  The  results  of  his  calculations 
agree  well  with  observed  soil  temperatures .  In  Experiment  1,  we  spec¬ 
ify  a  (Fig.  7)  as  a  function  of  surface  moisture,  represented  by  the 
ground  wetness  parameter  GW  in  our  model:  X  varies  from  1.2  *  10 
for  dry  soil  (GW  =  0)  to  6.0  ■  10  for  GW  >  0.2.  These  values  cor¬ 
respond  closely  to  those  given  by  Priestley  (1959).  The  depth  D  at 
which  the  diurnal  wave  of  soil  temperature  becomes  negligible  varies 
with  the  thermal  properties  of  the  soil.  In  a  dry  soil  the  tempera¬ 
ture  wave  does  not  penetrate  too  far,  whereas  in  a  moderately  wet  soil 
it  goes  much  deeper.  In  a  thoroughly  wet  soil,  D  is  the  one  nearest 
the  surface  (Chang,  1958).  Here  we  specify  D  to  vary  with  ground 
wetness  GW  so  that  D  increases  linearly  from  10  to  50  cm  as  the  ground 
wetness  GW  varies  from  0  to  0.5.  Then  D  decreases  linearly  from  50 


Ground  wetness 


Fig.  7 — Thermal  conductivity  as  a  function  of 
ground  wetness  used  in  Experiment  1 


to  10  cm  for  GW  =  0.5  to  1.  The  maximum  value  of  50  cm  is  chosen  on 
the  basis  of  observations,  which  show  that  for  common  soils  the  diurnal 
temperature  wave  does  not  penetrate  below  50  cm  (de  Vries,  1963)  .  We 


have  specified  TQ  to  be  280°K  at  all  bare  land  grid  points. 


Substituting  Eq.  (4)  in  Eq .  (3)  and  following  the  same  procedure 
as  in  the  control  experiment,  we  obtain  an  expression  for  computing 


ground  temperature  T  which  is  similar  to  Eq.  (2)  of  the  control  ex- 

O 


periment,  except  for  the  additional  terms  (X/D)T  in  the  numerator 


.  LJ 

and  /D  in  the  denominator,  which  arise  due  to  inclusion  of  the  soil 
heat  flux  in  the  heat  balance  equation  (3).  Therefore 
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Figure  8  shows  maximum  surface  temperature  distribution  for  the 
African  continent.  Comparison  of  this  with  corresponding  regions  in 
Fig.  2  (control  experiment)  shows  a  general  decrease  in  the  magnitude 
of  ground  temperature.  The  highest  value  of  surface  temperature  is 
7U°C,  as  compared  to  90°C  in  the  control  run.  The  maximum  amplitude 
of  the  diurnal  range  (Fig.  9)  also  shows  a  decrease  of  20°C  to  the 
west  and  10°C  to  the  east.  Other  regions  such  as  South  America  and 
Australia  also  show  an  improvement  over  control  experiment  results. 


Longitude 


Daily  maximum  ground  temperature  (°C) 
for  month  of  January  (Africa) 
(Experiment  1 ) 


Ground  temperature  (°C  )  Ground  temperature  (°C  )  Ground  temperatu 
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be  attributed  to  very  dry  surface  conditions  (GW  0)  at  locations 
A  and  B  but  completely  wet  conditions  (GW  0.9)  at  location  C.  Table 
3  shows  the  values  of  heat  balance  components  LE,  H,  and  G  relative 
to  the  net  radiation  at  locations  A,  B,  and  C  for  both  the  control 
experiment  and  Experiment  1. 


Table  3 


RATIO  OF  LE,  H,  AND  G  TO  NET  RADIATION  R 
AT  LOCATIONS  A,  B,  AND  C  N 


Locations  A 

,  &  B 

Location 

C 

Experiment 

LE/RN  1  H/RN 

G/RN 

1  LE/RN 

H/RN 

G/RN 

Control 

~  0  i  1 

| 

0 

1 

~  0 

0 

1 

~  0  0.6 

1 

0.4 

0.9 

~  0 

- J 

0.1 

tor  the  control  experiment  (G  =  0) ,  the  ground  temperature  at 
locations  A  and  B  (LE  0)  was  obtained  as  a  result  of  balance  betweer 
net  radiation  and  sensible  heat  flux  H,  whereas  at  location  C 
(H  0)  it  was  due  to  balance  between  R^  and  latent  heat  flux  LE. 

For  Experiment  1  (G  *  0) ,  at  locations  A  and  B  (LE  0)  G  is  about 
40  percent  of  net  radiation,  and  this  contributes  to  lower  tempera¬ 
tures  obtained  at  these  locations.  However,  at  location  C  the  balance 
is  still  predominantly  between  evaporation  and  net  radiation;  G  is 
only  10  percent  of  R^.  Consequently  there  is  no  significant  different 
between  ground  temperatures  obtained  at  C  for  the  control  experiment 
and  Experiment  1.  This  shows  that  while  heat  flux  into  the  ground  is 
important  for  dry  conditions  (  40  percent  of  Rn  for  bare  ground),  it 

is  not  very  significant  (  10  percent  of  R^  in  affecting  ground  tem¬ 

peratures  for  very  wet  conditions.  These  results  compare  well  with 
observations  (Sellers,  1965).  Just  as  in  the  case  of  the  control  run, 
Experiment  1  also  does  not  show  a  time  lag  between  maximum  solar  radia¬ 
tion  and  maximum  ground  temperature  (Fig.  11).  Similarly,  it  may  be 
noted  that  there  is  no  time  lag  between  maximum  soil  heat  flux  and 
maximum  ground  temperature.  These  results  also  follow  from  analytic 
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treatment  of  the  classic  heat  conduction  equation  if  we  neglect  heat 
capacity.  The  result  of  Experiment  1  is  at  variance  with  conditions 
in  nature,  which  show  that  the  maximum  surface  temperature  occurs  3 
hours  after  a  maximum  in  soil  heat  flux  (Johnson,  1954). 

Thus  Experiment  1  shows  that  while  the  inclusion  of  soil  heat 
flux  in  the  heat  balance  equation  gives  relatively  lower  values  of 
ground  surface  temperature,  the  relationship  of  ground  temperature 
with  solar  radiation  and  soil  heat  flux  is  not  realistic.  The  most 
important  drawback  of  the  formulation  of  G  used  in  this  experiment  is 
the  specification  of  one  value  for  1D  at  all  the  bare  land  points  on 
the  globe.  Moreover,  the  damping  depth  D.  which  can  vary  considerably 
over  the  globe,  also  has  to  be  specified. 


Experiment  2 

Gadd  and  Keers  (1970)  have  expressed  the  heat  flux  (G)  into  the 
ground  as  a  fraction  of  net  radiation  (B..)  at  the  surface  so  that 


=  U  (6> 


where  Rj,  -  (S  -  R) ;  and  u  is  assigned  a  value  of  0.1  for  L  j  0  (day¬ 
time)  and  0.5  for  RjJ  <  0  (nighttime).  Thus,  for  Experiment  2  we  ob¬ 
tain  the  ground  temperature  as  a  solution  of  the  equation 

H  +  LE  -  (1  -  u)rn  -  o  (7) 


(This  is  obtained  by  substituting  G  ■  uR^  in  Eq.  (3).) 

An  expression  for  computing  the  ground  temperature  T  (correspond¬ 
ing  to  Eq.  (5)  of  Experiment  1)  is  thus  given  by  8 


„t+l 


a  -  +  Vi  +  j„t  + 

r  ~T  <*q  a  )  ■ 
H1  ♦  r  G“  -jtH 


Tg  -  vv 


]( 


(8) 
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The  results  (fig.  12)  indicate  a  rather  insignificant  decrease  in  the 
daytime  magnitude  of  surface  temperature  as  compared  to  that  of  the 
control  experiment.  However,  during  nighttime  the  surface  temperatures 
of  Experiment  2  were  considerably  higher  than  those  of  the  control  ex¬ 
periment.  These  features  can  be  attributed  to  two  different  values  of 
“  used  in  Eq.  (2);  lower  for  daytime  and  higher  for  nighttime.  The 
values  of  diurnal  range  were  also  affected  only  to  the  extent  that  the 
minimum  in  Experiment  2  was  higher  than  that  in  the  control  experiment. 
It  may  be  mentioned  here  that  Gadd  and  Keers  (1970)  did  not  compute  LE 
and  H  explicitly,  as  we  did  in  Experiment  2.  They  obtained  the  latent 
and  sensible  heat  fluxes  by  partitioning  the  available  flux  (1^  -  G)  at 
the  surface. 

Thus  Experiment  2  does  not  improve  the  ground  surface  temperature 
distribution  obtained  in  the  control  experiment.  This  may  be  attributed 
to  the  technique  of  expressing  G  in  terms  of  R^p  which  gives  only  a 
crude  and  indirect  estimate. 

Experiment  3 

In  their  general  circulation  model,  Kasahara  and  Washington  (1971) 
computed  the  ground  surface  temperature  over  land  or  ice/snow  areas  by 
solving  Eq.  (3),  On  the  basis  of  a  study  by  Sasamori  (1970),  they  pre¬ 
scribed  the  soil  heat  flux  G  to  be  a  fraction  of  sensible  flux  H.  In 
addition,  they  specified  a  value  of  unity  to  the  Bowen  ratio  (H/LE). 
Thus,  if  we  substitute 


LE  =  H  j 


Eq.  (3)  becomes 


R  -  S  +  2.3H  =  0 


(10) 


or  substituting  for  H  in  terms  of  ground  temperature  T  ,  we  obtain 
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Fig.  12  Comparison  of  daily  ground  surface  temperature  (°C) 
between  control  and  Experiment  2 


-30- 
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The  results  show  that  the  values  of  the  ground  surface  temperature  in 
Experiment  3  are  much  smaller  than  those  obtained  in  the  control  ex¬ 
periment.  At  some  grid  points,  the  decrease  of  surface  temperature  is 
as  much  as  20  C.  Figure  13  shows  a  comparison  of  ground  temperature 
(at  two  locations)  between  Experiment  3  and  the  control  experiment. 

It  can  be  seen  that  while  the  maximum  has  decreased  considerably  (20°C) , 
the  minimum  has  increased  appreciably.  This  is  reflected  in  smaller 
values  (not  shown)  for  diurnal  oscillations  obtained  in  Experiment  3. 

Thus  Experiment  3  gives  a  much  improved  distribution  of  ground 
surface  temperature  vis-a-vis  the  control  experiment.  However,  the 
approximations  for  G  and  LE  (Eq.  (9))  are  not  quite  realistic.  For 
example,  the  first  condition  (G  =  0.3H)  implies  that  the  two  processes 
(G  and  H)  a»e  exactly  in  phase  and  have  the  same  constant  ratio.  This 
is  at  variance  with  observations.  Moreover,  the  value  of  0.3  for  this 
ratio  is  not  universal;  in  fact  Sasamori  estimated  this  value  for  dry 
soil  only.  Again  the  assumption  of  a  Bowen  ratio  of  unity  is  also  not 
realistic,  because  observations  (Brooks  and  Goddard,  1966)  have  shown 
that  it  has  a  significant  diurnal  variation  and  sometimes  changes  sign 
also.  It  may  be  mentioned  that  Kasahara  and  Washington  (1971)  suggest 
a  calibration  of  nighttime  surface  temperatures  by  controlling  the 
ground  temperature. 
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Yj _ PROGNOSTIC  equations  for  ground  temperature 

CALCULATIONS :  EXPERIMENTS  4  AND  5 

In  Sec.  IV  we  have  described  three  methods  of  computing  the  ground 
surface  temperature.  These  essentially  involve  solving  the  heat  balance 
equation  at  the  earth-atmosphere  interface  under  the  assumption  of  zero 
heat  capacity.  In  this  section  we  consider  the  situation  when  the  heat 
capacity  is  not  zero.  Table  4  shows  the  equations  used  for  solving 
f°r  T  in  Experiments  4  and  5. 

O 


Table  4 

EQUATIONS  USED  IN  EXPERIMENTS  4  AND  5 


Experiment 

Basic  Equation  to  Solve 
for  Ground  Temperature 

Expression  for  G 

4 

S-R-H-LE-G  =  0 

G  (0)  = 

(Tt+1  -  T  )  +  (Tt+1  -  TC  ) 

\  2 

\  8  /  tuAt  \  g  g/ 

3T 

r  •"  .  1 

\  9T  3 

5 

S-R-H-LE-G  =  c 

G(l)  = 

Hence  the  basic  difference  between  Experiments  4  and  5  is  that  the  for¬ 
mer  uses  the  heat  balance  equation  to  compute  T  ,  whereas  the  latter 
uses  an  explicit  predictive  equation  for  ground  temperature  T  .  Also, 
in  Experiment  4,  the  soil  heat  flux  G  is  evaluated  at  the  surface, 
whereas  in  Experiment  5  it  is  approximated  at  a  1-cm  depth  in  the  soil. 
Thus  while  the  equation  used  in  Experiment  4  involves  a  non-zero  heat 
capacity  only  implicitly  through  the  formulation  for  G;  the  equation  of 

Experiment  5  explicitly  involves  both  soil  heat  capacity  and  soil  heat 
flux  G. 

Experiment  4 

In  this  experiment  we  use  Eq.  (A-7)  of  the  Appendix  to  represent 
soil  heat  flux  (G)  at  the  surface.  Thus 
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Using  an  implicit  technique  and  writing  this  in  the  notation  used 
in  Eq.  (2)  of  the  basic  circulation  model,  we  have 


Here  u  is  the  frequency  of  oscillation,  c  is  the  volumetric  heat  ca¬ 
pacity,  X  is  the  thermal  conductivity,  T  is  the  average  daily  soil 
temperature  (assumed  to  be  the  s  une  at  all  depths),  and  At  is  the  time 
step  of  integration.  Substituting  this  expression  for  G  in  the  heat 
balance  equation  (3)  and  rearranging,  we  obtain 


3T 

- &  + 

to  Jt 


(12) 


(14) 


It  can  be  seen  that  Eq.  (14)  is  similar  to  Eq.  (2)  of  the  control  experi¬ 
ment  except  for  the  last  terms  in  both  the  numerator  and  denominator. 
These  arise,  in  a  way  similar  to  Experiment  1,  because  soil  heat  flux 
is  included  in  the  heat  balance  equation.  However,  for  the  control  ex¬ 
periment  and  Experiment  1  we  had  assumed  the  heat  capacity  c  to  be  zero; 
it  is  non-zero  for  Experiment  4.  For  X  we  assign  the  same  values  which 
were  used  in  Experiment  1  (Fig.  7).  As  regards  c,  it  can  be  expressed 
as  the  sum  of  heat  capacities  of  different  soil  constituents  in  a  unit 
volume  (de  Vries,  1963);  c  also  depends  upon  the  moisture  content.  In 
our  experiment  we  use  a  constant  value  of  0.6  cal  cm"3  deg"1,  which  is 
based  on  the  values  given  by  Priestley  (1959)  and  Geiger  (1965).  The 
frequency  of  oscillation  U)  is  7.27  x  10“5  sec"1,  and  At  is  30  min. 

Since  the  solution  for  primarily  varies  around  T,  the  accuracy  of 
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the  former  depends  to  a  great  extent  on  how  we  treat  the  latter.  The 
dally  average  T  can  be  Included  in  Eq.  (34)  In  two  ways: 

1.  Compute  T  on  the  basis  of  the  "past"  day's  values  of  T  and 
use  this  average  in  computing  Tg  for  the  following  day.  This 
technique  enables  us  to  let  T  vary  at  least  on  a  daily  basis; 
also  it  can  be  computed  at  each  grid  point. 

2.  Prescribe  T  as  a  function  of  latitude  zones  on  the  basis  of 
observations.  However,  since  there  are  no  maps  of  soil  tem¬ 
perature  available  on  a  daily  basis,  only  monthly  maps  can 
be  used.  Consequently,  solutions  for  Tg  for  any  month  on  a 

daily  basis  are  biased  toward  the  value  of  T  prescribed  for 
that  month. 


We  integrated  the  model  for  48  hours,  using  both  the  techniques  noted 
above.  We  did  not  find  any  significant  differences  in  the  results,  at 
least  for  the  two  days  we  considered.  The  results  discussed  here  are 
based  on  T  which  was  prescribed  as  a  function  of  latitude  zones  on  the 
basis  of  a  map  for  January  prepared  by  Chang  (1958)  (Fig.  14). 

Figure  15  shows  the  distribution  of  maximum  surface  temperature 
for  Experiment  4.  Comparing  this  with  the  results  of  the  control  ex¬ 
periment  (Fig.  2)  reveals  that  there  is  considerable  reduction  in  the 
magnitude  of  surface  temperature;  the  most  noteworthy  occurring  in  the 
western  half  of  Africa.  The  highest  value  in  this  case  is  only  60°C, 
wuereas  it  was  90°C  in  the  control  run.  Even  in  the  Sahara  region  the 
range  of  maximum  temperature  is  from  30e  to  50eC;  for  the  control  run 
these  values  are  from  30°  to  70°C  with  50e  and  60e  isotherms  dominating 
the  region.  In  Australia,  South  Asia,  and  South  American  continents 
also,  the  values  tend  to  become  more  realistic.  It  may  be  noted  that 
in  this  case,  as  also  in  Experiment  1,  there  are  still  rather  unreal¬ 
istically  high  values  of  surface  temperature  in  the  equatorial  forest 
region  of  Africa.  The  reason  for  this  is  that  during  the  duration  of 
the  integration  of  the  model  (48  hours)  there  has  been  no  precipitation 
in  this  region  and,  quite  contrary  to  a  "real"  atmospheric  situation, 
there  is  no  moisture  to  cause  evaporr tional  cooling.  Simulation  of 
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January  climate  by  Cates  (1972)  has  shown  this  region  to  be  character¬ 
ized  by  heavy  precipitation  and  consequently  a  high  evaporation  rate. 

It  can  perhaps  be  stated  that  this  unrealistic  feature  of  relatively 
higher  surface  temperatures  in  the  equatorial  African  region  will  dis¬ 
appear  in  any  long-term  integration  of  the  model. 

Figure  16  shows  distribution  of  the  diurnal  range  of  ground  temp¬ 
erature.  Again  the  unrealistic  features  displayed  by  the  control  run 
(Fig  3)  have  improved  considerably,  especially  in  Australia,  South 
Asia,  South  America,  and  Africa.  Figure  17  shows  a  plot  of  T  ,  G,  and 
solar  radiation  S  as  a  function  of  local  time.  It  can  be  seen  that 
there  is  a  phase  lag  of  about  3  hours  between  maximum  G  and  maximum  T  , 

g  ’ 

and  a  lag  of  about  1  hour  between  maximum  S  and  maximum  T  .  These 

g 

features  agree  very  well  with  observations  (Sellers,  1965;  Lettau, 
1951),  as  well  as  analytical  solutions  of  the  heat  conduction  equation, 
and  are  produced  by  the  interaction  of  various  components  of  the  energy 
balance  at  the  surface. 

Experiment  5 

Recently  some  atmospheric  circulation  models  have  computed  the 
ground  surface  temperature  (T  )  by  solving  an  explicit  prediction  equa¬ 
tion  for  Tg.  According  to  Corby  et  al.  (1972),  the  ground  surface 
temperature  can  be  computed  from  the  equation 

3T 

C  3^  *  S  "  R  "  H  "  LE  (15) 

where  c  is  the  heat  capacity  per  unit  area  of  the  surface  layer,  and 
where  S,  R,  F  and  LE  have  the  same  meaning  as  in  the  control  and  other 
experiments.  It  may  be  noted  that  Eq.  (15)  does  not  incorporate  soil 
heat  flux  G  explicitly.  Consequently,  we  are  not  able  to  reproduce  the 
diurnal  variation  of  the  surface  temperature  (T  )  realistically.  In 

this  experiment  we  describe  our  attempt  to  incorporate  soil  heat  flux  G 
in  Eq.  (15). 

As  shown  by  Sellers  (1965),  for  example,  in  the  absence  of  hori¬ 
zontal  temperature  gradients,  the  rate  at  which  a  soil  layer  is  warming 
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Diurnal  range  of  ground  surface  temperature 
(°C  )  -  Experiment  4  (January  ) 
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or  cooling  is  proportional  to  the  net  flux  of  energy  into  or  out  of 

the  layer.  Thus  for  a  layer  of  thickness  Az  the  observed  time  rate 
of  temperature  change  is  given  by 


-c 


AT 

7^  Az  =  AG 


or,  in  differential  form 


c 


^T 
_ s 

at 


aG 

3z 


(16) 


Also,  as  shown  in  the  Appendix,  the  soil  heat  flux  G  can  be  expressed 
as  a  function  of  z  and  t.  Thus 


G(z.t) 


(oct 


+  COS  (u)t 


-i>] 


(17) 


or 


G(z,t)  *= 


U  3T 


+  Ts(*,t) 


(18) 


In  practice  it  is  preferable  to  use  Eq.  (18)  rather  than  Eq.  (17) 
because  ATq,  the  amplitude  of  temperature  wave  at  the  soil  surface,  Is 
rarely  known  (Sellers,  1965).  De  Vries  (1963)  and  others  have  related 
ATq  to  the  amplitude  of  air  temperature  at  2  meters  for  the  annual 
cycle,  but  such  a  relationship  cannot  be  determined  for  the  daily  cycle 
because  of  the  great  variability  involved. 

Let  us  consider  a  soil  layer  extending  from  the  surface  (z  *  0)  to 
a  depth  z.  Then  by  Eq.  (16), 

3T  (z) 

— §£ - -  lG(z)  -  G(0) ]  (19) 
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where  Tg(z)  is  the  average  soil  temperature  of  layer  of  depth  z 

G(0)  =  S-R-H-LE  Is  the  soil  heat  flux  at  the  surface  z  =  0  and 
Involves  ground  temperature  (T  )  of  R,  H,  and  LE,  and 
dependence 

G(z),  (Eq.  (18)),  Is  the  soil  heat  flux  at  depch  2. 

However,  since  we  were  Interested,  mainly,  1„  computing  the  ground 
surface  temperature  Tg,  and  also,  since  our  model  had  no  explicit  levels 

within  the  soli,  we  used  Eqs.  (18)  and  (19)  to  obtain  an  equation  that 

enabled  us  to  predict  T  . 

g 

Let  us  consider  a  soil  layer  of  1-cm  thickness.  Since  there  Is 
considerable  uncertainty  In  defining  "surface,"  particularly  on  global 
scale,  and  also  because  the  horizontal  grid  size  of  our  model  is  quite 
crude  (4°  latitude,  5°  longitude),  we  can  assume  that  average  soil  tem¬ 
perature  for  a  layer  of  1-cm  depth  approximates  the  ground  surface 
temperature  itself.  Thus 


T  (1)  i  T 
s  •  g 

Applying  this  approximation  to  Eqs.  (18)  and  (19)  at  z  -  1  cm,  we 
obtain 


G(l)  * 


f] 


and 


9T 

=  -  IG(1)  -  G(Q)  ]  =  G(0)  -  G  (1) 


Thus  we  have 


ST 

= 

at 


S-R-LE-H- 
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3T 


u.  a  t 


+  T 
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or,  rearranging, 


3T 


'1  at 


S  -  R  -  LE 


(T  -  T) 
§ 


(20) 


where 


c 


1 


It  may  be  noted  that,  except  for  S,  all  terms  on  the  right-hand 

side  of  Eq.  (20)  are  functions  of  ground  surface  temperature  T  .  Thus 

g 


=  F(T  )  =  S  -  R 
g 


LE  -  H 


"  T) 


Using  the  Newton-Rapbston  technique  and  appropriate  formulations 
for  R,  LE,  and  H  in  F(T  )  above, 

O 


T(t+1)  = 
g 


T 

g 
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S  -  R  -  H  -  LE 


~  +  4vT3  +  C..U  + 

-t  g 


'-ui1 


(21) 


Here  =  c  +  .  c/2~  (c  being  volumetric  heat  capacity,  prescribed  as 

0.6) 

=  thermal  conductivity  (same  as  in  Experiment  1) 
r  =  average  daily  surface  temperature,  assumed  same  at  all  depths 
as  prescribed  i..  Experiment  4 
*  frequency  of  oscillation  (same  as  in  Experiment  4) 
t  =  time  step  of  integration  =  30  minutes  (1/48  day) 

=  Stefan-Boltzmann  constant  =  1.171  •  lo”7  ly  day”1  deg”^. 


Other  terms  have  the  same  meaning  as  in  the  earlier  experiments  of  Sec. 
IV.  The  terms  on  the  right-hand  side  are  to  be  evaluated  at  time  t. 
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We  used  Eq.  (21)  to  predict  ground  surface  temperature  for  a  48- 
hour  period  of  January.  Figures  18  and  19  show  global  distribution  of 
daily  maximum  ground  surface  temperature  and  its  diurnal  range  for 
January.  A  comparison  of  this  with  the  control  experiment  shows  a  sub¬ 
stantial  decrease  in  the  magnitude  of  the  ground  temperature.  The  re¬ 
sults  of  Experiment  5  look  similar  to  those  of  Experiment  4  (Fig.  15) 
except  for  the  African  continent,  where  the  regions  of  high  maximum 
surface  temperature  (>60'C)  are  much  smaller  in  Experiment  5  than  in 
Experiment  4.  As  regards  diurnal  range  (Fig.  19),  the  amplitudes  are 
smaller  by  10°  than  those  shown  by  Experiment  4  (Fig.  16)  in  the  equa¬ 
torial  African  region.  The  phase  relationships  between  T  ,  S,  and  G 
(not  shown)  are  the  same  as  those  demonstrated  ir.  Experiment  4:  the 
maximum  in  Tg  follows  the  maximum  in  solar  radiation  (S)  by  about  an 
hour  and  follows,  by  about  3  hours,  the  maximum  in  soil  heat  flux  G. 

Even  though  results  of  Experiments  4  and  5  appear  to  be  similar, 
it  may  be  noted  that  Experiment  5  uses  an  explicit  prediction  equation 
for  Tg,  whereas  in  Experiment  4  it  is  implicit  in  the  formulation  of 
G.  It  may  also  be  remarked  that  whereas  the  formulation  of  Experiment 
4  is  applicable  only  at  the  surface,  that  of  Experiment  5  is  more  gen¬ 
eral  and  can  be  used  to  predict  the  temperature  of  a  soil  layer  of 
finite  depth  and  consisting  of  multiple  levels. 


LONGITUDE 


maximum  ground  surface  temperature 
Experiment  5  (January) 


LONGITUDE 


Diurnal  range  of  ground  surface  temperatu 
°C)  -  Experiment  5  (January) 
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VI.  SUMMARY  AND  CONCLUDING  REMARKS 


In  accordance  with  the  purpose  of  this  study,  we  have  integrated 
the  Rand  two-level  general  circulation  model  to  compute  the  ground 
(bare  land)  surface  temperature  (Tg)  by  incorporating  soil  heat  flux 
in  the  heat  balance  equation  (Experiments  1  through  3).  We  have  also 
computed  Tg  by  taking  into  consideration  the  heat  capacity  of  the  soil  as 
well  as  soil  heat  flux  (Experiments  4  and  5).  A  comparison  of  results 
(Table  5)  shows  that  the  most  realistic  distribution  of  T  with  respect 
to  the  magnitude,  the  diurnal  range,  and  the  phase  relationships  is  ob¬ 
tained  in  Experiments  4  and  5.  However,  since  the  formulation  for  Experi¬ 
ment  5  is  more  general  and  can  be  used  to  compute  the  surface  temperature 
of  a  soil  layer  of  finite  depth,  it  should  be  preferred  to  that  of  Experi¬ 
ment  4,  especially  if  one  can  afford  a  model  of  multiple  soil  levels. 

The  technique  of  Experiment  4  is  considered  adequate  if  only  the  inter¬ 
face  temperature  is  required.  Other  methods  such  as  the  control  ex¬ 
periment,  Experiment  2  (Cadd  and  Keers,  1970),  and  Experiment  3  (Kasa- 
haru  and  Washington)  are  not  realistic;  the  latter  two  involve  extensive 
parameterization  of  the  soil  heat  flux.  This  study  has  also  demonstrated 


Table  5 


COMPARISON  OF  THE  RESULTS  OF  VARIOUS  NUMERICAL  EXPERIMENTS 
AT  LOCATION  B  (IN  SAHARA  REGION) 


Ground  Temperature  (°C) 

Experiment 

Maximum 

Diurnal 

Phase  Relationship  Between  Maximum 

Range 

Values  of  Tg,  S,  and  G 

Control 

40 

46 

1 

34 

35 

(  Maxima  in  Tg,  S,  G  occur  at  same  local 
(  time 

2 

39 

38 

3 

32 

31 

) 

4 

31 

28 

1  Maxima  in  Tg  occur  1  hour  after  maxima 

5 

—  i 

27 

-  1 

22 

_ 1 

in  S  and  3  hours  after  maxima  in  G 

NOTE:  Tg  is  ground  temperature,  S  is  solar  radiation,  G  is  soil  heat  flux. 


a  good  correspondence  between  the  results  of  numerical  experiments  and 
the  classic  analytical  solutions. 

We  have  not  considered  the  temperature  of  surfaces  other  than  bare 
land.  In  order  to  simulate  natural  conditions  we  should  also  consider 
surfaces  which  may  be  covered  with  ice  and/or  snow.  However,  the 
atmospheric  model  used  in  this  study  prescribes  a  cover  of  ice  and/or 
snow  which  does  not  change  with  time.  Therefore  we  have  not  modified 
Eq-  (1)  for  computing  ground  temperatures  of  surfaces  other  than  bare 
land.  Modifications  of  these  conditions  are  currently  in  progress. 
Besides  the  nature  of  the  cover  on  the  surface,  we  should  also  consider 
its  other  properties  such  as  texture,  moisture  content,  organic  matter, 
color,  structure,  cultivation,  and  topography. 

The  results  described  in  this  study  are  based  on  only  a  48-hour 
integration  of  the  basic  circulation  model.  Consequently  the  conclu¬ 
sions  cannot  be  generalized  for  long-term  integration  results. 

Since  we  have  integrated  a  basic  circulation  model  for  different 
experiments  after  varying  only  the  formulations  for  soil  heat  flux, 
this  study  presents  a  valid  comparison  of  the  usefulness  of  different 
methods.  It  is  felt  that  the  techniques  discussed  can  also  be  used  in 
integrating  local  atmospheric  and  soil  boundary  layer  models. 

In  recent  times,  considerable  attention  has  been  paid  to  the  study 
of  climatic  changes  and  their  impact  on  worldwide  economic  and  social 
order.  Though  modeling  of  climatic  variation  is  still  in  developing 
stages,  general  circulation  models  are  being  used  increasingly  to  study 
shorter  period  variations.  And  recognizing  the  importance  of  feedback 
effects  of  variations  of  surface  temperature  on  atmospheric  pro^sses 
it  is  essential  to  be  able  to  compute,  among  other  things,  a  realistic 
distribution  of  ground  surface  temperature.  The  technique  recommended 
in  this  report  can  be  useful  in  this  computation. 
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Appendix 

AN  EXPRESSION  FOR  HEAT  FLUX  INTO  THE  SOIL 


Assuming  that  soil  is  homogeneous  and  that  heat  flow  is  in  the 
vertical  direction  only,  the  heat  conduction  equation  can  be  written  as 


3T  ,  32T 

_ s  A^  _ s 

3t  c  2 

3  z 


(A-l) 


where  T&  is  soil  temperature,  A  is  thermal  conductivity,  and  c  is  volu¬ 
metric  heat  capacity. 

Let  us  assume  that  at  all  depths  z  the  temperature  varies  as  a  pure 
harmonic  function  around  an  average  value.  Though  crude,  this  assump¬ 
tion  enables  us  to  approximate  the  description  of  actual  fluctuations 
caused  by  the  succession  of  day  and  night  (diurnal)  or  winter  and  sum¬ 
mer  (seasonal).  However,  this  must  be  modified  for  natural  soil  temper¬ 
atures.  Let  us  now  assume  that  the  temperature  at  the  surface  can  be 
written  as 


T  (0, t)  =  T  +  AT  sin  (wt)  (A-2) 

o  O 

where  T  is  daily  average  temperature  of  the  soil  assumed  to  be  the  same 
at  all  depths,  ATq  is  the  amplitude  at  the  surface,  and  U)  is  frequency 
of  oscillation  and  equals  2A/ (period  of  the  wave).  Equation  (A-2)  is 
the  boundary  condition  at  z  =  0.  The  solution  of  (A-l)  may  be  written 
as 


Tg(z,t)  =  T  +  ATQe"z/d  [sin  (wt  -  z/d) ]  (A-3) 

where  T  (z,t)  is  the  soil  temperature  at  the  depth  z  and  time  t  and 
d  =  «/2A/cuj  is  the  depth  at  which  the  amplitude  of  ATq  is  insignificant. 

Again,  for  an  infinitesimally  thin  soil  layer,  the  her.t  flux  into 
the  soil  is  given  by 


Preceding  page  blank 
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9T 

G(z,t)  =  -X  ■—  (A-4) 

Combining  (A-3)  and  (A-4),  we  get 

G(z,t)  «=  ATo  e  Z//d  [sin  (wt  -  z/d)  +  cos  (cot  -  z/d)]  (A-5) 

Also,  from  Eq.  (A-3)  we  obtain 

ATQe  z^d  sin  (cot  -  z/d)  * 
and  differentiating, 

ATQe  Z ^  cos  (cot  -  z/d) 

Therefore,  by  eliminating  ATq,  which  is  rarely  known  accurately,  Eq. 
(A-5)  can  be  written  as 


Tg(z,t)  -  T 


co  3t 


G(z,t) 


/SfilVlIi 
V  2  u  3t  + 

P 


Tc(«.o 


(A-6) 


or  at  z  =  0  (surface) 


G (0 , t) 


(A-7) 


where  T  -  T  (0,t)  is  the  temperature  at  the  surface. 

O  ^ 
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